Conductance oscillations in metallic nanocontacts 
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We examine the conductance properties of a chain of Na atoms between two metallic leads in the 
limit of low bias. Resonant states corresponding to the conductance channel and the local charge 
neutrality condition cause conductance oscillations as a function of the number of atoms in the 
chain. Moreover, the geometrical shape of the contact leads influences the conductivity by giving 
rise to additional oscillations as a function of the lead opening angle. 
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I. INTRODUCTION 

The prospects to fabricate molecular devices, extend- 
ing from a few atoms to 10 /im long carbon nanotubes, 
have inspired an active research field Q . The most mod- 
ern experimental, theoretical and computational meth- 
ods are used to investigate and model ultimate-size tran- 
sistors, switches, gates, and memory devices as well as 
their coupling with the rest of the electronic circuit. 
In this work we concentrate to chains of Na atoms 
as ultimate conductors. Similar atomic-scale metallic 
nanowires have been produced from Au by a scanning 
tunneling microscope (STM) §§) and by the mechani- 
cally controllable break junction (MCBJ) technique jjj. 
The formation of Na nanowires has been investigated by 
first-principles calculations ||-f7j and the ensuing conduc- 
tances have been compared with experiments [[r| . 

The conductance of a nanowire shows nonmonotonic 
behavior as a function of its length. However, the role of 
the contact leads cannot be ignored, as shown by Sim et 
al. H and Yeyati et al. ||. The purpose of this article 
is to study the dependence of the conductance not only 
on the length but also on the shape of the contacts. We 
show that the shape of the contact, i.e. the bluntness of 
the STM tip, influences the conductance in an oscillatory 
way. 

Yeyati et al. f| recognized the importance of resonance 
states to the conductance of sharp leads. These reso- 
nances exist only at sharp tips, not at blunt ones. In the 
case of monovalent atoms, such as Na or Au, the appear- 
ance of a resonance at the Fermi level can explain the 
first quantized conductance step as a multiple of 2e 2 jh. 
Sim et al. M showed that the charge neutrality required 
for the Na chain and the tips of the contacts leads to an 
even-odd behavior in the conductance as a function of 
the number of atoms in the chain. This is because for 
an odd number of atoms the resonance is half-occupied 
whereas for an even number of atoms the resonance is 
fully-occupied, leading to a low density of states at the 
Fermi level. 

As Sim et al. || we see the even-odd behavior for lead 
tips. In addition, the conductance oscillates as a function 
of the lead opening angle. This results from the change 
in the angular character of electrons screening the tip 



regions of the leads. 



II. MODEL 



As in Ref j8j , we model the nanocontact by a chain of 
N Na pseudoatoms |n| . To enable a systematic study of 
the shape effects, we model the leads using the "stabilized 
jcllium" model |H| ], where leads of discrete ions have been 
smeared out to continuous positive charge distributions 
and a constant potential energy term has been added 
to stabilize the Na bulk electron density. The computa- 
tional geometry is depicted in Fig. [j]. The positive back- 
ground charge of each cone corresponds to M electrons. 
The constant density n + of the positive background is 
determined from the Na valence electron density param- 
eter r s — 3.93 a.u. (n + — 3/(47rrf)). Near the atom 
chain, the radius of the cones is r s and the larger radii 
are determined by the number of electrons in the cones, 
which can be varied at will. We solve the electronic states 
of the system using the density-functional theory within 
the local density approximation |12| assuming the limit 
of zero external bias. 




FIG. 1. Geometry of the present model. Na pseudoatoms 
are located between two jellium cones. The cone angle a can 
be varied continuously. The dotted line gives the boundaries 
of the calculation cell. For the details, see also the text. 

The system has cylindrical symmetry, so that the prob- 
lem is computationally two-dimensional. In the middle 
of the atom chain there is a mirror plane perpendicular 
to the symmetry axis. The electronic wave functions are 
required to be even or odd with respect to the mirror 
plane. Naturally, the Poisson equation involved in the 
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self-consistent loop is solved by requiring mirror symme- 
try for the Coulomb potential. On the other boundaries, 
i.e. on a plane parallel to the mirror plane but far outside 
the jellium cone and on a cylinder surface far from the 
axis, the wave functions are required to vanish. The nu- 
merical solutions of the density-functional equations are 
obtained by a real-space multigrid scheme |l3[ . 

In the case of a chain of Na atoms there is only one con- 
ducting eigenchannel contributing to the electron trans- 
port at the Fermi level. Then the mirror symmetry and 
the Friedel sum rule give the conductance G in the zero- 
bias limit as |l3] 



G = —— sui 2 
h 



N ) 



(1) 



where N e and N Q are the numbers of electrons in the even 
and odd parity states, respectively. In an infinite system, 
N e and N Q are calculated by integrating the density of 
states up to the Fermi level. The states are assumed to 
have the spin degeneracy and to be occupied accordingly. 
Then the difference of N e and iV arises only from the 
electrons of the conducting eigenchannel, because only 
their wave functions may have a finite value at the mir- 
ror plane. In our model the conducting eigenchannel cor- 
responds to the quantum number to = of the angular 
momentum around the symmetry axis. 

The states with a finite value at the mirror plane 
(to = 0) span an energy region starting somewhat be- 
low the Fermi level. They are strongly localized at the 
chain atoms and the lead tips. These resonance states, 
as occupied up to the Fermi level, have to participate in 
neutralizing the positive charge of the atom chain and 
the tips of the leads. According to the calculations by 
Sim et al. ||, for Na atom chains between two Na (111) 
tips, the neutrality requirement leads for an odd number 
N of Na atoms to a half-filled resonance in the conduc- 
tion eigenchannel. This means that the Fermi level is at 
the resonance state, which is filled in a spin-compensated 
manner totally by one electron. As a result, the difference 
N e — N Q is odd and the conductance is close to 2e 2 /h. 
For an even N the Fermi level is off the resonance and 
the conductance will be less than unity. Below we will 
show that also electron states which do not correspond 
to the conducting eigenchannel (to > 0) participate in 
the charge neutralization in the tips, modifying the con- 
ductance behaviour. 

In order to study the convergence of the results as a 
function of the system size we have performed calcula- 
tions with M = 200 and 300 electrons in each jellium 
cone. The finite size results in a discrete eigenvalue spec- 
trum which is smoothed using the Fermi-Dirac broaden- 
ing with a finite temperature when occupying the elec- 
tron states. The broadening should be larger than the 
single electron level spacing but smaller than the distance 
between the resonances. The level spacing decreases as 
~ M- 2 / 3 so that we can decrease the temperature with 
increasing M and obtain sharper resonance structures. 



For 200 electrons we have chosen the temperature of 600 
K, which was used also by Sim et al. 0] for cones of 60 
and 95 electrons. In the case of 300 electrons in a cone 
we choose the temperature of 400 K. 



III. RESULTS 

Our main result is the finding that the conductance of 
a chain of atoms depends on the cone angle of the leads. 
In order to understand the reasons we first discuss the 
changes in the electronic structure in terms of such in- 
tegrated quantities as the electron density and the local 
density of states (LDOS). Thereafter the conductance, 
which reflects changes in the electronic structure on a 
microscopic level, is given. 




FIG. 2. Electron densities induced by chains of one, two, 
three and four pseudoatoms between two jellium leads. The 
cone angle is 69° and each of the cones contains 300 electrons. 
The density values of the contours are given in atomic units 
(a~ 3 ). Areas of positive induced density are shown as gray 
and the ionic positions as black circles. 

Fig. |^ shows the electron density induced by chains 
of one to four Na atoms between two jellium leads con- 
taining 300 electrons each. The cone angle is 69°. The 
induced density is calculated as the difference between 
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the whole system and the corresponding system consist- 
ing of the leads only. In accordance with the calculations 
M having 95 pseudoatoms as leads, the main part of the 
induced density is localized at the atom chain and the 
density decays rapidly in the cones. This indicates con- 
vergence as a function of the system size. For an even TV 
the electron density is slightly more localized than in the 
case of an odd N because for an even N strong resonance 
states are below the Fermi level and are fully occupied. 
The density of states (DOS) of the system with two Na 
atoms is shown in Fig. ^. It is obtained from the discrete 
energy level spectrum given at the bottom of the figure 
by substituting the lines by Lorentzians, the widths of 
which correspond to the temperature used in the self- 
consistent electronic structure calculations. We see that 
around the Fermi level the single electron level spacing 
is small in comparison with distances between the strong 
DOS peaks. The latter are related to the m = res- 
onance structure (the m — eigenenergies are given in 
the figure as long vertical lines). It can be seen that the 
region of peaks with a larger amplitude reaches about 1 
eV below the Fermi level. 
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FIG. 3. Chain of two Na atoms between two jellium leads 
with the cone angles of 69°. The density of states (DOS) is 
given with respect to the Fermi level (energy zero). The single 
electron eigenenergies are marked with bars at the bottom of 
the figure. The m = states are denoted by longer bars than 
the m > states. 

The character of the electron density near the tips of 
the leads changes as the cone angle increases. This is 
demonstrated in Fig. || which shows the total electron 
density and its m = state component for three dif- 
ferent cone angles. The densities are integrated over a 
plane perpendicular to the cylinder axis and are given as 
a function of the position z along the axis. Naturally, 
the total electron density increases close to the tip as the 
cone angle increases. But it is important to note that 
the m > electron states contribute most of the density 
increase. As a result, the m = resonance states feel an 
increased Coulomb repulsion and are pushed upwards in 



energy. For the same reason, also other states with low 
m (> 0) values rise in energy, but to a lesser extent. 
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FIG. 4. Chain of two Na atoms between two jellium leads. 
The electron density integrated over the plane perpendicular 
to the cylinder axis as a function of the position z along the 
axis is given. The upper panel shows the total density whereas 
the lower panel gives the contribution due to the m = states, 
only. Densities for the three cone angles of 42°, 60°, and 80° 
are given by solid, dotted and dashed lines, respectively. The 
vertical line marks the jellium edge and the black circles de- 
note the Na atom positions. 

The shift of the m = states to higher energies can 
be seen in the local density of states (LDOS) at the Na 
atom chain. We calculate the LDOS similarly to the 
DOS discussed above, but now the area of a Lorenzian 
peak is the density of the state in question integrated 
over a volume limited by two planes perpendicular to the 
symmetry axis. In the left panel of Fig. || we show the 
LDOS of the even and odd m — states for the cone an- 
gles of 67°, 68°, and 69°. The integration is limited to the 
atom chain, i.e. between the planes at the central jellium 
edges. The LDOS of the chain has a structure similar to 
the result of the atomistic lead calculations [|| . When we 
open the cone angle, the m — states are shifted, as dis- 
cussed above, to higher energies. The conductance of the 
chain is close to 2e 2 /h when the Fermi level is located on 
top of a peak of the m = LDOS. The left panel of Fig. 
^| indicates that the conductance varies rapidly as a func- 
tion of the cone angle. The right panel of Fig. || gives 
the total LDOS calculated by including the tips of the 
cones into the integration; the limiting planes are at the 
depth of r s from the central jellium edges. We see that in 
this volume, which is important for the charge neutrality 
of the atom chain, the m = features of the left panel 
are hardly visible. The change of this total LDOS is also 
less dramatic than that of the m = LDOS at the atom 
chain. 
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FIG. 5. Chain of two Na atoms between two jellium leads. 
The local density of states (LDOS) near the Fermi level (en- 
ergy zero) is given for different cone angles. The left panel 
shows the the even (solid line) and odd (dashed line) m = 
LDOS's calculated for the atom chain between the jellium 
edges. The right panel shows the total LDOS's calculated for 
the atom chain and the tips of the leads to the depth of r 3 
from the jellium edge. The LDOS's corresponding to the two 
uppermost angles are shifted in steps of 100 and 200 units in 
the left and right panels, respectively. 

On the basis of Figs. ^ and [| we can conclude that 
the to > states have an important contribution to the 
local charge neutrality and influence the conductance of 
the chain. This is seen in Fig. ^| which gives the conduc- 
tances of the chains of one to four Na atoms correspond- 
ing to these three cone angles. The even-odd oscillation 
as a function of the number of chain atoms is clear as 
well as its gradual phase change. 
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FIG. 6. Conductances of the chains of one, two, three and 
four Na atoms between two jellium leads. The conductances 
are shown for three different jellium cone angles. The number 
of electrons in each cone is 300. 

Finally, Fig. |?] shows the conductance for the N = 3 



chain as a function of the cone angle. Results for system 
sizes corresponding to 200 and 300 electrons in each jel- 
lium cone are given. It can be seen that for both system 
sizes the conductance is a regularly oscillating function 
of the cone angle. The temperature is closely connected 
to the amplitude of the oscillations so that a lower tem- 
perature gives stronger oscillations. If we choose a high 
enough temperature, the oscillation as function of angle 
disappears but so does the oscillation as function of N. 
The number of oscillations seems to increase as the sys- 
tem size increases. This is because the resonance struc- 
ture becomes more resolved including higher to compo- 
nents as the size of the system increases. 
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FIG. 7. Chain of three Na atoms between two jellium 
leads. The conductance is given as function of cone angle for 
systems of 200, and 300 electrons in one lead. The tempera- 
ture used in occupying the one-electron states is 600 K for the 
smaller system and 400 K for the larger system. The dotted 
line shows the conductance of the chain of two atoms. 

The differences in the curves for the systems of 200 and 
300 electrons in Fig. (7] indicate that the finite system size 
has still some influence. However, there are also clear 
similarities. For example, there appears shallow conduc- 
tance minima between deeper minima. This structure is 
connected to the clustering of the resonance m = en- 
ergy levels as seen in the level spectrum at the bottom 
of Fig. |. The effect can be considered as reminiscent of 
the shell structure in the system of conical confinement. 

In Fig. ^ we show for the larger system the conduc- 
tance of the chain of two Na atoms at large cone angles. 
The conductances of the N = 2 and N = 3 chains are 
in opposite phases and according to our calculations this 
state continues at least to the angle of 110°. Thus, we 
were not able to find a limiting blunt angle at which the 
influence of the resonance states and the even-odd oscil- 
lations in the conductance disappear . 

The oscillation of the conductance as a function of 
the cone angle means the alternation of perfect electron 
transport and partial electron scattering in the to = 
channel. The electron scattering is caused by the interfer- 
ence between the to = channel and the non-conducting 
to > channel. The situation is analogous to the trans- 
port through a quantum wire with a side-coupled quan- 
tum dot |ll],[l(| . If the Kondo resonance of the quantum 
dot lies at the Fermi level the transport channel DOS at 
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the Fermi level is suppressed (the situation is also called 
as the Fano antiresonance) and the conductance has a 
dip. In the present model, the m > channels play the 
role of the side-coupled quantum dot. 

As discussed in the Introduction, atomic chains have 
been experimentally observed for Au . These chains 
have been interpreted to contain up to four Au atoms. 
It has been noted that before breaking the last conduc- 
tance plateau of about one conductance quantum is not 
smooth but exhibits abrupt changes synchronized with 
abrupt changes in the elongation force JlTj]. The steps 
in the conductance may be of the order of one tenth of 
the conductance quantum. Our findings give two expla- 
nations for these small conductance changes. The con- 
ductance may jump abruptly when the number of atoms 
in the chain increases. On the other hand, the geometry 
of the lead tips changes affecting also the conductance. 
These two effects may be superimposed during the actual 
elongation process. 

IV. CONCLUSIONS 

We have studied the electronic properties of chains 
of Na atoms between two metallic leads using a 
pscudoatom-jcllium model. The conductances of the 
chains are calculated using the Friedel sum rule. The con- 
ductances show even-odd oscillation as function of num- 
ber of atoms in the chain confirming the earlier results 
obtained by a purely atomistic model There is an- 
other oscillation: the conductance oscillates as a function 
of the lead cone angle. The even-odd oscillation is a direct 
result of the resonance structure of the electron states of 
the conducting channel and the charge neutrality at the 
atom chain and the tips of the leads. The conductance 
oscillations as a function of the cone angle reflect the fact 
that also electrons of the non-conducting channel take 
part in the charge neutralization and that the balance be- 
tween the conducting and non-conducting channel elec- 
trons depends on the cone angle. Thus, when determin- 
ing the conductance of a nanoconstriction it is important 
to solve self-consistently for the electronic structure, in- 
cluding the lead regions. 
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